Frequency-dependent admittance of a short superconducting weak link 
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We consider the linear and non-linear electromagnetic responses of a nanowire connecting two 
bulk superconductors. Andreev states appearing at a finite phase bias substantially affect the 
finite-frequency admittance of such a wire junction. Electron transitions involving Andreev levels 
are easily saturated, leading to the nonlinear effects in photon absorption for the sub-gap photon 
energies. We evaluate the complex admittance analytically at arbitrary frequency and arbitrary, 
possibly non-equilibrium, occupation of Andreev levels. Special care is given to the limits of a single- 
channel contact and a disordered metallic weak link. We also evaluate the quasi-static fluctuations 
of admittance induced by fluctuations of the occupation factors of Andreev levels. In view of possible 
qubit applications, we compare properties of a weak link with those of a tunnel Josephson junction. 
Compared to the latter, a weak link has smaller low-frequency dissipation. However, because of the 
deeper Andreev levels, the low-temperature quasi-static fluctuations of the inductance of a weak 
link are exponentially larger than of a tunnel junction. These fluctuations limit the applicability of 
nanowire junctions in superconducting qubits. 
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The search for longer coherence times of supercon- 
ducting qubits brought the study of finite-frequency elec- 
tromagnetic properties of mesoscopic superconductors to 
the forefront of experimental research— ~—. The majority 
of experiments until recently was performed on structures 
using Josephson junctions as "weak" superconductors, 
and substantial progress in recognizing the coherence- 
limiting mechanisms was achieved. One may view a num- 
ber of mechanisms causing energy or phase relaxation as 
extrinsic ones. These involve, e.g., imperfections in the 
tunnel barriers comprising junctions 7 -, charge trapping^, 
and interaction with stray photons*^. Along with them, 
there are intrinsic mechanisms associated with the kinet- 
ics of quasiparticles in the superconductors^ - — . These 
mechanisms provide fundamental limitations to the co- 
herence. The majority of effects of quasiparticles on the 
finite-frequency properties of Josephson junctions can be 
derive d 12 ' 14 from the electromagnetic admittance of the 
junction Y(uj). This property was extensively studied 
theoretically, starting from the seminal phenomenologi- 
cal paper of Josephson— and microscopic evaluatio n 16 ' 17 
based on the BCS theory. 

The use of weak superconducting links instead of 
Josephson junctions in qubits was proposed recently as 
a way to avoid extrinsic decoherence mechanisms (such 
as imperfections of the tunnel barriers)^. An apparent 
observation of a coherent phase slip in a conducting weak 
linki^ may be viewed as an incipient experimental step in 
that direction. That makes the question about the intrin- 
sic mechanisms of decoherence in weak links important. 
Like with Josephson junction devices^, this question is 
directly related to the finite-frequency admittance of a 
weak link. Surprisingly, this property was given relatively 
little attention to. The admittance of a short SNS con- 
tact was investigated, mostly numerically, in the recent 
paper o 20 ' 21 . Some qualitative aspects of the AC response 
of a single-channel point contact can be extracted from 
two other papers devoted to the theory of enhancement 
of supercurrent by microwave radiatio n 22 ' 23 . 



Here we perform a fully-analytical evaluation of the 
admittance of a weak link connecting two bulk supercon- 
ductors, valid at arbitrary frequency oj, quasiparticle dis- 
tribution function, and normal-state conductance of the 
link. Compared to the Josephson junction case, the dis- 
sipative part of the weak link admittance exhibits a num- 
ber of new thresholds in its frequency dependence, asso- 
ciated with the presence of Andreev levels. The complex 
admittance close to these new threshold frequencies is 
sensitive to the occupation of the discrete Andreev states. 
Fluctuations of the equilibrium or non-equilibrium occu- 
pation factors result in fluctuations of the admittance. 
We analyze the average values and fluctuations of the 
linear electromagnetic response, giving special attention 
to the practically important limits of a single-channel 
contact 2 ^ and a disordered metallic wire 3 -^. 

The discrete nature of Andreev states is responsible 
for a low threshold for the nonlinear absorption. In the 
nonlinear regime, we find a suppression of the absorption 
coefficient in a disordered metallic link at radiation fre- 
quency ui < 2 A/ 3, while at higher u dissipation power 
depends non-linearly on the radiation intensity (here A 
is the BCS gap in the leads). 

The paper is organized as follows: the model used in 
the derivation of the admittance of a point contact with 
an arbitrary transmission coefficient is formulated in Sec- 
tion HI The linear response theory for the AC perturba- 
tion of the point contact is developed in Section [TTJ In 
Section Mil we discuss the results for the admittance of 
the point contact at zero temperature and no quasipar- 
ticles present. In Section IIVI we study the changes in 
the admittance caused by the arbitrary distribution of 
quasiparticles in the junction. These results are used in 
Section [V] to find the admittance of a disordered weak 
link. The fluctuations of the admittance are analyzed in 
Section IVI1 both for the case of point contact and of a 
weak link. In Section rVIII we consider the absorption rate 
in a non-linear regime for the radiation frequencies close 
to the Andreev level resonance. We conclude with the 



2 



final remarks in Section [Villi 



I. POINT CONTACT HAMILTONIAN 

We start by considering a point contact between two 
leads. It can be described by the tunnel Hamiltonian 



H — Hl + Hr + Ht ■ 



(1.1) 



in V requires at least the smallness of the perturbation 
to the dynamics of the system, |</>i| = |eV/cj| -C 1, where 
u> is the frequency of perturbation. Further limitations 
on the parameters, which may come from the effect of V 
on occupation factors, will be discussed later. 

It is convenient to do the gauge transformation Ck — > 
Cfce 1 ^ 1 before performing the perturbation theory. This 
moves the </>i-dcpcndcnce to the tunneling terms. Using 
the Kubo formula for linear response, we get 



where -Hl(.r) are the BCS Hamiltonians of the left (right) 
leads: 

Hl = J2 k &4<* + A l E fe 44 + A l E fe c- k c k , (1.2) 
H R = J2 P £p4 c p + A R S P 44 + A *R c p c p , (1-3) 



and 

H T = (4 c p + 4 c ^) + W J2 4i c *a + w 

kp kik 2 P1P2 

(1.4) 

is the tunneling term. Here c&(p) and ^ are elec- 
tron operators in the left (right) lead corresponding to 
states k(p) and its time-reversed pairs k{p), and Al(R) = 
Ae'^w) are the BCS gap functions. 

The tunneling amplitude w is assumed to be 
momentum-independent near the Fermi level. It is re- 
lated to the transmission coefficient, r = {2-kvqw) 2 /[I + 
(2m/Qw) 2 ], where vq is normal-state density of states. 
The conductance of the junction G in the normal state 
is proportional to r (hereinafter we set K = 1). A point 
contact between superconducting leads hosts a single An- 
dreev level with energy Ea{t, <f) depending on G: 

G = e 2 r/7r, Ea(t, <j>) = A(l — r sin 2 <f>/2) 1 ^ 2 ; (1.5) 

here the phase difference between the leads order param- 
eters, <j) — <j)R — (j>L, is assumed to be time-independent. 



(i{t)) = i,,+ dt' x (t-t')Mt') 



(2.3) 



Here Ij is the Josephson current which is present even 
without applied voltage: 



I j = ew Im I 4^p 

kp 

The response function x(t) is given by 



(2.4) 



X (t) = iew 2 8(t) £ J2 { [4, (*)dpi (*) ~ < (*) c *i W< 

fclPl fc 2 p 2 

4 2 (0)d P2 (0) - 4 2 (0)c fe2 (0)] ) - ewS(t) Re ( ]T 4d p ) . 

kp 

(2.5) 

Averages (...) are taken over the Gibbs ensemble of the 
original Hamiltonian H . We can use Wick's theorem to 
evaluate averages. They can be expressed in terms of 
Green's functions of the unperturbed system. Green's 
functions satisfy a system of linear integral equations, 
but the corresponding kernels are separable due to the 
form of the tunneling term Eq. (| 1 .4[) . Therefore, that 
system reduces to a system of algebraic equations which 
can be solved exactly. This response function is related 
to the admittance in frequency domain: 



Y(u>) = i-xM ■ 



(2.6) 



II. LINEAR RESPONSE TO AC 
PERTURBATION 

We may account for an applied small, time-dependent 
voltage V(t) by modifying <Pl -> 4>l + 2(f>i(t) in Eq. (|1.2[) , 
with 0i = eV(t), and adding the term — eV(t)N~L to 
Eq. dTT]): 



ft = H(t) - eVN L 



Nr 



E4< 



We want to find the current (J), 
/ = eN L 



/cp 



(2.1) 



(2.2) 



induced by an applied voltage to linear order in V and at 
arbitrary transmission r. The validity of linear response 



It is convenient to split the admittance into a sum 



Y 



Y 



loLj 



(2.7) 



each term of which has a clear physical origin. The purely 
inductive term Lj comes from the u — > response of the 
condensate, 



1 

Lj 



2e 



01. 



J 



(2.8) 



The other five parts of Eq. (|2.7|) originate from the quasi- 
particles transitions. To better understand the structure 
of these parts, recall that ReY (dissipative part of ad- 
mittance) is related to the linear absorption rate W of 
the radiation by: 



W 



2e 2 



uj\ReY(uj) 



(2.9) 



3 



The elementary processes leading to the absorption are 
depicted in Fig. 1. The Yi(uj) term corresponds to a pro- 
cess in which two quasiparticles are created in the band, 
leading to the energy threshold 2 A. The contribution 
Y2{uj) in (|2.7|l comes from creating one quasiparticle in 
the bound state and one in the band, the correspond- 
ing threshold energy is A + Ea- Creation of a pair of 
quasiparticles in the bound stated, which costs energy 
2Ea, leads to the term Y3 (u>). In addition to these three 
contributions which exist even in the absence of quasipar- 
ticles, there are two more associated with the promotion 
of an existing quasiparticle to a higher energy in the ab- 
sorption process: Y4 is the intra-band contribution, and 
Y5 corresponds to an ionization of an occupied Andreev 
level. 



III. ADMITTANCE OF A SINGLE-CHANNEL 
JUNCTION AT T = (EQUILIBRIUM STATE 
WITH NO QUASIPARTICLES) 

Evaluating averages in Eq. (|2.4|) for system at zero tem- 
perature we get for the Josephson current: 



(0) 



A 2 sin § 
2eE A {4>) 



(3.1) 



in agreement with the result one obtains from Eq. ()1.5[) 



by differentiating energy over <p, i.e. / 



(0) 



-2edE A /d(j). 



Using the above expression for Ij and Eq. (|2.8jl we find 
the inductive term at zero temperature: 



(0) 



ttGA- 



COS ( 



^ sin 4 I 



^ sin 2 ± 



3/2 



Evaluating averages in Eq. (|2.5[) we find the contri- 
butions Y^ to the admittance. With no quasiparticles 
present, there can be no processes of type 3 or 4 in Fig- 



ure [TJ Therefore, y} 0) (lo) = Y^ v) (uj) = 0. The contribu- 



te), 



tion K L (0) (w) comes from the creation of pairs of quasipar- 
ticles in the band (two excitations of type 1 in Figure [TJ, 
and its real part is given by: 



RcY^V) _0(w-2A) 
G ~~ u 



;-A 



dep(e)p(tL> - e)\z e 



(3-3) 



where p(u>) is the density of states in the continuum nor- 
malized to normal-state density of states z/q: 



eVe 2 - A 2 
e 2 - El 



(3.4) 



We assume u > throughout this section. The result 
for negative frequencies can be found using the fact that 
ReY^°'(uj) is an even function. The theta function in 
Eq. (|3.3p shows that there can be no creation of pairs in 
the continuum for frequencies less than 2 A. 

The Y 2 (oj) term comes from processes in which one 
quasiparticle is created in the band and another one in 
the Andreev level. These processes are represented by 
one arrow of type 1 and one of type 2 in the Figure [TJ 
The real part of F 2 (w) is given by: 



Rey 2 (0) M 

G 



= tt9{lu-E a -~A) 



-E 2 

—p{lo-E a )\z Ea ^\ 2 



(3.6) 



and it vanishes for u < A + Ea, as for these frequencies 
the processes "1+2" are energetically not allowed. 

Finally, there are processes in which two quasiparti- 
cles on Andreev level are created. Those are represented 
by two excitations of type 2 in Figure [TJ In this case, 
the frequency must be equal to 2Ea- The Y^°\uj) term 
comes from such processes and its real part is given by: 



Rcy 3 (0) H 

G 



(A 2 -£ 2 )(£ 2 -A 2 cos 2 f) 



2E\ 



l ' S(u)-2E A ) . 

(3-7) 



Note that the RHS of Eqs. (j33|)-([3Jl) depend on G and 
cj) through Ea, see Eq. (|1.5j) . 

The admittance exhibits non-analytical behavior at 
threshold frequencies ui = 2Ea, A + Ea and 2 A. For 
u> w 2A we have ReF 1 (0) oc (w - 2A) 2 (9(w - 2 A) ac- 
(3-2) cording to (|3.3j) . Similarly, for frequencies u> ~ A + Ea 
we get ReF 2 (0) cx y/w — (A + E A )9(u - (A + E A )) from 
Eq. (EU). 

The imaginary parts of li(a;)'s can be obtained from 
their real parts using Kramers-Kronig relations since 
Y(lu) is analytic in the upper half of the complex cj-plane. 
The complete expression for Imy(w) is given in Ap- 
pendix A. At threshold frequencies Im Y(w) exhibits non- 
analytical behavior which parallels threshold behavior of 
Re Y(w). Atw~ 2 A the non-analytical contribution be- 
haves as ImVj ' cx (lu — 2A) 2 ln|w — 2A| and at uj ps 

A + Ea it behaves as ImF 2 (0) cx y/(A + E A ) - uj9((A + 
Ea) — w)). The coefficients omitted from the asymptotes 
of Re7^(w) and lmY^ \oj) equal each other, confirm- 
ing that the complex function Y^(uj) is analytical. 



IV. ADMITTANCE OF A SINGLE-CHANNEL 
JUNCTION IN THE PRESENCE OF 
QUASIPARTICLES 



and the matrix element z is given by: 



ja = 1 _ A»ooe* + A»-lft (3 _ 5) 



e(w - e) 



The admittance changes once there are quasiparticles 
present. Each term in Eq. (|2.7p acquires an additional 
factor depending on the quasiparticle occupation num- 
bers. We introduce occupation factors po, pf, p± and p 2 
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FIG. 1: Schematic representation of contributions to the ad- 
mittance, Eqs. (|3.3[) - (|3.7p . (|4.5l) and <|4- 7 P - Horizontal bars are 
energy levels and arrows depict possible excitations. The pic- 
ture shows that Ya = Ys = in the absence of quasiparticles. 
Since ReY(w) is proportional to the absorption rate, we see 
that Re Vi, 2,3,5 will have threshold frequencies of 2A,A + _E.4, 
2Ea, A — Ea, respectively. 



denoting probabilities of having zero, one or two quasi- 
particles in the bound state; po + + pi + P2 = 1. The 
inductance in Eq. (|2.7[) then becomes: 



1 

Lj 



L 



1 

To) 



bo - P2] 



(4.1) 



The Y\ term acquires a factor depending on the occupa- 
tion factors /(e) of the continuum states, 



Relf^u;) _0(w-2A) 



G 



u—A 



dep(e)p(uj - e)\z e ^\' 



x [1 - /(e) - /( W - e)] . 



(4.2) 



This expression is different from Eq. p.3[) by a factor- 
equal to the difference of probabilities for having the ini- 
tial and the final state occupied in the transition from 
the ground state to the band, see Fig. Q] Similarly, the 
Re Yi term is given by 

Re Y 3 (u) = Re Y 2 °M(Po + - f(u> - E A )) , (4.3) 



and Y3 term by 



Y 3 (u) = Yg(u)(p - Pa ). 



(4.4) 



At non-zero occupancies, there are two additional con- 
tributions to absorption, Y^{u) and Y^w). The for- 
mer one comes from the band-to-band transitions, repre- 
sented by the arrow 4 in Figure [T] Rs real part is given 
by: 



RcYijcu) _2 
G ~w 



dep(e)p{u + e)\z.^\ 2 [f{e)-f{uj + e)] 



The other term, Ys(w), is generated by the Andreev level- 
to-band transitions^ 2 -. These transitions are represented 



by the arrow 3 in Figure [TJ We can write it in the form 
resembling that of Y^fa): 



ReY 5 (iu) = ReF 5 (0) H[p t + Pi + 2p 2 - 2f(u+E A )} 



(4.6) 



where the real part of Y^ (u) is given by: 



Rer 5 (0) H 
G 



'-tt6{uj+E a -A) 



2lo 



■p(uj + E a )\z_ Ea ^\ 2 , 
(4.7) 



The occupation factors p and / in all of the above ex- 
pressions may, but need not to be the equilibrium ones. 
The density of states p(e) and matrix clement z c . w in 
Eqs. (|4~2"j) and (|4~5j) - P~7)l are defined in Eqs. (pT4"]) and 
(EH). 

Note that the transitions involving the Andreev level 
vanish at E A — > A. This is achieved if either G = or 
= 0. Expanding to the lowest order in G for G <C e 2 /7r, 
the contributions Y1.4 are linear in G and they reduce to 
the familiar perturbative result for the admittance of a 
Josephson junction 2 ^. The other contributions are higher 
order in G, with Y 2 , 5 oc G 3/2 , and ReF 3 cx G 2 . 

In the limit 0—^0, the terms involving transitions 
to Andreev level vanish as Re Y^s oc \<fi\, and ReY"3 cx 
<f> 2 . Therefore, at tf> = only terms contributing to Y 
are again Y\ and Y4. In that case, the expression for 
Y coincides with the one found perturbativelj*21 in the 
limit of weak tunneling G <C e 2 /7r. Thus, at small phase 
bias </> we don't expect much difference from the simple 
Josephson junction. 

Now we analyze the behavior of admittance in the lim- 
its of low frequencies and low temperatures, as these are 
the conditions often encountered in the application of 
the superconducting junctions. At frequencies below the 
threshold for Andreev level ionization, u < A — E A , and 
away from the bound pair creation resonance, u> 7^ 2E A , 
the only contribution to the dissipative part of admit- 
tance comes from Eq. (|4.5[) . We assume the quasiparti- 
cle occupation factors are distributed according to Boltz- 
mann distribution /(e) = Ae~ e / T . At low temperatures 
T <C A — E A , where A — E A is characteristic scale for 
the energy dependence of the density of states above the 
gap, the dominant contribution to Re Y^ioS) comes from 
the transitions between the states near the bottom of the 
band. In that limit, we get for the asymptotic form of 
ReY 4 (o;): 



ReY" 4 (w) 




(4.8) 



(4.5) where x qp is the density of quasiparticles, n qp , in the bulk 



normalized to the "Cooper pair density" , x 



''ii' 



and Ki (x) is the modified Bcsscl function of the second 
kind. Note that at small frequencies, u> <C T, it follows 
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from Eq. (|4.8[) that Re 1-4(0;) is frequency- independent 
and proportional to y/T/ A. 

In the limit G <C e 2 /7r, the Andreev level is shallow, 
A — Ea -C A. If now T S> A — £U, the main contribution 
to Y4 comes from transitions involving states far above 
the gap where the density of states is described by the 
usual BCS result. In this limit, Eq. (|4.5[) is reduced to 
the known result^ for a Josephson junction, 



Rer(w) 
G 



1 



2 % /2^oA 



(1 



COS( 



3/2 



AT 



= In — • 

II u 



(4.9) 



This limit is the opposite to the one of Eq. (j4.8[) . The 
two asymptotes match each other at T ~ A — Ea up to 
the logarithmic factor. 

It is interesting to compare the dissipation in a large- 
area Josephson junction of G < e 2 /tt with the dissipation 
in a single- channel weak link of the same G. The weak- 
link quasiparticle density of states in the continuum, see 
Eq. (|3.4p . is suppressed compared to the singular tun- 
neling density of states in a Josephson junction. As a 
result, at frequencies and temperatures u>,T < A — Ea 
a weak link is less dissipative than a Josephson junction 
with small-transparency but large-area tunnel barrier of 
the same G. Using Eqs. (|4.8[) and (|4.9p we find, e.g., 
that the dissipation is smaller by a factor T 2 /A 2 in the 
case of a weak link. 



V. DISORDERED WEAK LINK 

For a multi-channel junction one needs to sum the con- 
tributions to the admittance from each channel. We 
consider the case of a disordered weak link for which 



we can assume the transmission coefficients arc contin- 
uously distributed according to Dorokhov distribution 2 ^ 
p(r) = -KG/2e 2 T\Jl — t. The admittance is then given 
by: 



drp{T)Y(u,T) . 



(5.1) 



We can write Y(u>) as a sum of five terms, in the same 
way we did it for the single-channel junction in Eq (|2.7[) . 
Transitions between Andreev levels are ignored, which is 
justified in the limit of short junction with A/ Et — > 0, 
where Et is Thouless energy 2 ^. 

The Josephson current of the disordered weak link 
can be found using the same averaging procedure as in 
Eq. (|5.ip . In the absence of quasiparticles it is given by: 



-(o) ttGA <p . , 

i>= cos — artanh sm ■ 

J e 2 



Similarly, averaging 1/L { j ] from Eq. we get: 



:(0) 



= ttGA 



1 — sin — art anh sin — 
2 2 



(5.2) 



(5.3) 



Evaluating the integral in Eq. (|5.1[) results in expressions 
for Yi. If there are no quasiparticles present, the only 
non- vanishing terms are Y± 2,3- Their real parts exhibit 
threshold behavior at frequencies 2 A, A + A| cos ^| and 
2A|cos-||. The latter two thresholds correspond to the 
fully transmitting channel which has the lowest possible 
Andreev level energy at given phase <p- The term corre- 
sponding to the creation of a pair of quasiparticles in the 
continuum is given by: 



RegV) = _ 2A)p yT du)i "I" - A 2 (l + cos0) ^Jj^Ky^-Lo.Y- A 2 ^ y^/A 2 - cos 2 \ + | sin || 
G 7 1 A 2 W 2 |sin||(w-2w 1 ) ^ W 2/ A 2 _ C0S 2 1 ^f/A 2 - cos 2 | - | sin || 

(5.4) 



It has the threshold frequency 2 A, just like the single 
channel admittance, Eq. (|3.3|) . At frequencies higher 
than this threshold, ReYi(oj) starts to grow linearly, 



ReYi(w) aw - 2 A. The ReY^u;) term, correspond- 
ing to the creation of one quasiparticle in the Andreev 
level and one in the continuum, is given by: 



RelfV f . . ujloi — A (1 + cos <p) y^£^ ^A^y^ 2 ^ 2 

= tt / duiBiu - <Ji - A) : — — . (5.5) 

J w 2 A|smf|(w-2 Wl ) X 2 -A 2 cos 2 | 

A|cos|| V 1 2 



I 

The threshold frequency of this term is A + A| cos t\, same as the threshold frequency of the fully transmitting 



6 



channel for this process. Behavior near threshold is given 
by Re Y 2 (lu) oc to — A — A | cos ^ | . Finally, there is a term 
coming from the processes in which a pair of quasiparti- 



J 



cles is created in the Andreev level. It has the threshold 
frequency of u> t h (</>) and is given by: 



Rer 3 (0) H 

G 



= 0(w - w t/i (0))0(2A - wV 



A 



sin f| w 



1 - 



4A 5 



r 



(0) 



Uth{4>) = 2A|cos-| 



(5.6) 




FIG. 2: Real part of admittance of a weak link. Solid lines 
correspond to cj> = 7r case. The explicit formulas are given by 
Eqs. (|5.4[) . (|5.5|) and (|5.6j) . Re Y 2 and Re Yi exhibit threshold 
behavior at uu — A and w = 2A where they start to grow as 
oc (u~ A) 3 ^ 2 and oc lu — 2A, respectively. Re Y3(ci>) is diverging 
for — > 0. Dashed line shows Re Y?,(u>) at phase different than 
7r, with threshold frequency ui t h, see Eq. (|5.6|) . As — !> 7r its 
maximum grows and shifts towards the lo — 0. 



In the presence of quasiparticles, the above expressions 
for the admittance acquire additional factors reflecting 
the quasiparticles distribution function, similar to the 
single-channel junction. In addition, there are two other 
terms, Y4 and Y5, coming from band-to-band transitions 
and ionization of Andreev level, respectively. These are 
obtained by averaging Eqs. (|4.5j) and (|4.6[) over transmis- 
sion coefficients. The complete expression for the dissi- 
pative part of the admittance in the presence of quasi- 
particles is given in Appendix B. 

We expect the greatest change in admittance from the 
simple Josephson junction at <f> « 7T, when Andreev level 
energies of the channels contributing to the admittance 
fill the whole range of energies between and A. In that 
case, for uj t h < lu < A and no quasiparticles present, the 
only contribution to dissipative part of the admittance 
comes from Y^{lu) and is given by Eq. (|5.6p with Lu t h — 
A\4>— 7r|. At this threshold Re Y-^w) grows as [ui — uith) 1 ' 2 
and reaches maximum for u> = y/2u>th- The height of 
the maximum scales as A/uj t h- Frequency dependence of 
RcY 3 for <p close to 7r is shown on Fig. 2. 



When <j) = 7T exactly, there is no low-frequency cut- 
off. For low frequencies ReY^w) diverges as l/u>. At 
frequencies A < oj < 2 A, in addition to RCY3, there is 
the contribution Rel^ given by Eq. (|5.5[) . Its behavior 
near threshold for <j) = tt is different than for any other <p 
and is given as Re Y2 oc (w — A) 3 / 2 . At frequencies higher 
than 2 A, the ReYSj contribution vanishes and ReY = 
ReYi -l-Rel^. The frequency dependence of ReY(oj) for 
cf> = n is also shown on Fig. 2. 

If </> ^ 7r, the dissipative part of the admittance is zero 
for tu < tuth and vanishing occupation factors. Assum- 
ing Boltzmann distribution f(E) = Ae~ E l T for quasi- 
particles and considering frequencies lu < LUth, the only 
contribution to ReY(u;) comes from I4 and I5 terms 
- Eqs. (|B6|) and (|B7I) . The former comes from transi- 
tions within the continuum band. In the limit lu,T <C 
A(f — I cos I ) the most important are transitions from 
the bottom of the band and we get: 



G 



2tt 



cot 



%V<%), (5.7) 



where U(x) = dte~ xt ^t{l + t)ln(l + l/t). The 
Re Yc,{lu) term is due to transitions from Andreev levels 
to the continuum. In the same limit of small frequency 
and temperature, it is given by: 



Rcy 5 (u; 



G 




r AT , to lu , . 

M^f), (5-8) 



where Ii(x) is the modified Bessel function of the first 
kind. The dissipative part of admittance is given by the 
sum of the two terms: ReF(w) = ReY^w) + ReY^w). 
For lu -C T the leading term comes from Eq. (|5.7|) and is 
frequency- independent, Re Y4 oc \J A/T. Comparing the 
considered case of a weak link to a tunnel junction of the 
same conductance G, ReY^w) has an additional, small 
factor of T/(Asin 2 § ln(T/w)), suggesting that Re Y 4 (lu) 
is reduced in the case of a weak link at low frequencies. In 
the opposite limit, lu T, the leading term comes from 
Eq. (|5.8[) due to higher population of low energy Andreev 
levels. In that case, Re Y 5 (lu) oc x qv e^l T ^TATlu^I' 2 . Be- 
cause of the large exponential factor, the dissipation in 
the weak link is greater than in the tunnel junction of 
the same conductance at high frequencies. 



7 



VI. FLUCTUATIONS OF ADMITTANCE 

Superconducting junctions are crucial elements of su- 
perconducting qubits. The admittance of a junction af- 
fects the properties of such qubits (i.e. their frequency 
and relaxation rates )i 2 -. As shown above, the admittance 
depends on the number of quasiparticles in the junction. 
Fluctuations of the occupation numbers cause fluctuation 
of the admittance. Consequently, the resonant frequency 
of a qubit containing the junction will fluctuate. 

The inductive 1/Lj term (which determines the fre- 
quency of the qubit) depends only on the occupation 
numbers of Andreev level. Therefore, its variance de- 
pends only on the variance of the occupation numbers of 
Andreev level 



(4°>) 



i Var ^- = (p + p 2 - (p - P2) 2 ) , 



.1 



also Eq. (|5.3|) . Under the same assumptions, the variance 
of the mean inductance 1/Lj is: 



L ( r 0) Vvar 1 



val Lj ~ 2irG Xqpe 
(0 = 



(i-|co S ||)A/T. ff((/))) (65) 



cos|| 5 / 2 



sin f |(1 — sin * artanhsin |) 2 



(6.6) 



The factor e 2 /2irG can be interpreted as 1/N e , where 
N e is the effective number of channels in the weak link. 
Comparing to the case of a weak tunneling junction with 
the same number of channels, the relative fluctuations of 
1/Lj have a factor exp((l — | cos ^|)A/T). This large ex- 
ponential factor suggests that the fluctuations are greater 
in the weak link. Such shot-to-shot fluctuations con- 
tribute to inhomogeneous broadening and limit the use- 
fulness of weak links in superconducting qubits. 



The relative fluctuations of Ij and 1/Lj are significant 
unless po , P2 or p^ + pi are close to 1 . Assuming equi- 
librium between quasiparticles in the band and in the 
Andreev level, as well as low quasiparticle occupation 
numbers of the Andreev level, Eq. (J6.1D reduces to: 



4°)y Var J_ = : 



qp 



(A/T) 



1/2(A-E A )/T 



(6.2) 



VII. NON-LINEAR ABSORPTION RATE AT 
RESONANT FREQUENCY 

For frequencies w ~ 2Ea we found that the admit- 
tance of a single-channel junction has a resonant delta- 
function peak corresponding to creation of quasiparticles 
at the Andreev level. Using Eq. (|3.7[) , we may re-cast the 
absorption rate Eq. (|2.9I) in the form 



The Yi terms in the expression (|2.7p depend on the oc- 
cupation numbers of the continuum states as well. How- 
ever, the fluctuations of admittance caused by the fluctu- 
ations of these occupation numbers are inversely propor- 
tional to the volume of the system and therefore are neg- 
ligible in the macroscopic limit. The expression for the 
variance of Y(lu) is then similar to the one for Var l/Lj. 
At frequencies uj < A + Ea and lu 7^ 2Ea to avoid the 
resonance, we get: 



VarRey(w) w (p + p 2 



bo-p 2 ) 2 )(RcF 5 (0) M) 2 , 

(6.3) 



where KeY^°\u;) is given by (|4.7p . Note that at low 
frequencies ReY 5 (uj) = as it has a phase-dependent 
threshold. 

To calculate the fluctuations in disordered weak links 
one must integrate the above expressions for variances 
over the distribution of the transmission coefficients. As- 
suming again equilibrium between band states and An- 
dreev levels, the fluctuations of the Josephson current are 
given by: 



Var^^AsJ^'^ 1 -'^ 
2 qp I ^ 41 



u e (l— |cos f |)A/T _ 



(6.4) 



Here we also assumed T <C A(l — | cos ^|), so that the 
main contribution comes from the channels with lowest 
Ea (fully transmittive channels), and T <C A| cos t\, see 



W = ^n 2 R 6{u) - 2E A ) . 



Here 



n, 



(A 2 — E\) 



E\ - A 2 cos 2 § 



AE 



a sin ■ 



(7.1) 



(7.2) 



has the meaning^ of Rabi frequency for the transitions 
in an effective two-level system driven by AC perturba- 
tion (fix. The two levels correspond, respectively, to the 
empty and doubly-occupied Andreev state. In the lin- 
ear response, we neglect the effect of Rabi oscillations on 
the dynamics of the two-level system. This is possible as 
long as fi ft is smaller than some "natural" , independent 
of 4>i width rjQ of the levels. Such natural width coming, 
e.g. , from inelastic scattering of quasiparticles 2 ^ leads to 
a replacement S(u — 2Ea) — > r] /iv[(u} — 2Ea) 2 + rj 2 ]. 

The effect of a stronger AC perturbation is two-fold. 
First, it may result in > rj affecting the dynamics 
of the two-level system. Second, it may make the levels 
lifetimes dependent on <f>\ by introducing new processes 
in the kinetics of quasiparticles. Indeed, the AC field may 
"ionize" the Andreev state, transferring a quasiparticle 
from that state into the continuum. One needs lu > A — 
Ea for that. Using the resonance condition, u 2Ea, 
we find that the kinetics of the Andreev states is sensitive 
to the AC perturbation at Ea > A/3. 

To address these two effects, we truncate the time- 
dependent part of the original Hamiltonian (|2.1[) retain- 
ing only terms responsible for the Rabi oscillations be- 
tween the empty and doubly-occupied Andreev state and 



8 



terms causing the ionization of that state, 

a ka 

+ a|aj] + cosu)t^\Xka\ a aq + X* k a\a k a] ■ (7.3) 

ka 

Here a a and a ka arc annihilation operators of quasiparti- 
cles in the Andreev level and the band, respectively. The 
last sum in Hamiltonian (|7.3|) is responsible for the tran- 
sitions between the Andreev state and continuum. The 
corresponding ionization rate is 



k \ 2 S{uj 



E A )- 



(7.4) 



The very same term leads to the Re Y5 part of admittance 
in the linear response theory, allowing us to relate r/ to 
Rey 5 , 



n 



^E A ReY 5 {0 \2E A ) . 



(7.5) 



The Hamiltonian (|7.3[) is quadratic, so the equations 
of motion for operators a reduce to a linear system of 
differential equations. Assuming frequencies close to the 
resonance, u — 2Ea « A - Ea, we can find the behavior 
of the solutions to the equations of motion after a long 
period of time, t 1/rj. The system then reaches the 
stationary state in which the energy absorption rate P is 
given by: 



P = (H) = -n A uri , 



(7.6) 



with tia being the average number of quasiparticles in 
the Andreev level in the process of Rabi oscillations, 



n A 



O 2 



(w - 2E A ) 2 + n 2 R + if 



(7.7) 



Hereinafter we neglected a shift of the resonant frequency, 
\2Ea — 2Ea\ oc (f>i, which is paramctrically smaller than 
the broadening due to the Q R + r] 2 term in the denomina- 
tor of Eq. (|7.7[) . The expression for the absorption power 
Eq. (|7.6[) has a simple interpretation: 77 is the transition 
rate from the level to the band, so uaij is the rate at 
which the Andreev level loses quasiparticles. To keep 
the number of quasiparticles in the level stationary, for 
each particle that left, a new one must be created in the 
level. This amounts to energy of uj + Ea = 3uj/2 for each 
transition, explaining the factor of 3/2 in Eq. (|7.6p . The 
condition Ea > A/3 needed for 77 7^ and the resonant 
condition oj ~ 2Ea imply uj > 2 A/3 for the absorption 
power Eq. (|7.6p to be finite. 

Using Eqs. (HTJ), (EH), g2D, iH, (& and (1731) 
for Haj ?7, and £Ir we can write P in terms of the Rabi 
frequency, 



37iE A n 2 R 



[OJ 



2E A f 



Mr 



V = 



n W 9E A - A2 ( E A + A 2 COS 2 |) 

16E^/W = E A \ 7 E~ 



- A 2 cos 2 f ) 



(7- 



At generic values of static phase bias <f> and transmission 
coefficient t, one has flu 3> 77 as long as the perturbation 
is reasonably weak, cf>i <C 2ir. In that case ua exhibits 
saturation at resonance, while P grows linearly with the 
perturbation intensity oc 4> 2 . At a small static phase bias 
77 takes form 



V8 



(7.9) 



It indicates that an increase of the excitation amplitude 
may result in a non-monotonic ha vs. ^Ir dependence 
and in saturation of P at fairly low excitation strength, 
J2 fl ~ |0| 3 A^(l~r). 

The population of the Andreev states by quasiparticles 
drastically alters the critical current of the junction and 
its low-frequency properties due to the changes in the 
inductance. Using Eqs. (|4.1I) and (|7.7|) we find 



1 



1 



(w - 2£ A ) 2 



Lj L f (uj-2E A ) 2 + n 2 R + V 2 



(7.10) 



Therefore, Qr may be inferred experimentally from a 
measurement of the critical current^! or from a two-tone 
experiment of the type^. 

In the case of a disordered weak link, an AC voltage at 
frequencies 2A| cos ^| < u> < 2 A also populates Andreev 
levels with quasiparticles. If the applied voltage is low 
and the ionization processes are negligible, levels with 
energies within an interval \Ea — oj/2\ < are substan- 
tially populated, cf. Eq. (|7.7|) . The resulting absorption 
power, 

3 40e^ 
x [iw 2 + A 2 cos <f> + A 2 



w 2 A 2 sin 2 I 



(7.11) 



scales as \4>i\ 3 reflecting the growing with |0i| number of 
states involved in the absorption. As before, it is required 
that uj > 2 A/3 to allow the AC- field- induced ionization 
of the excited Andreev levels. 



VIII. CONCLUSION 

The motivation for this study was two-fold. First, it 
came from the prospects^ of using nanowires instead 
of tunnel junctions in qubits and other microwave de- 
vices^. Additional impetus for the study came from ex- 
periments^ with nano-scale junctions, pointing to their 
extreme sensitivity to the presence of quasiparticles. 

We obtained an analytical expression for a frequency- 
dependent admittance of a point contact of arbitrary 
transmission coefficient and arbitrary quasiparticle oc- 
cupation factors. The results are valid even for non- 
equilibrium distribution of quasiparticles (see Section 
lIVp . The generalization to a short weak link (shorter 
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than coherence length), is presented in Section [V] We 
found that at low frequencies and temperatures, which 
are of interest in qubit devices, the dissipation of a point 
contact and a disordered weak link may indeed be lower 
than in a tunnel junction of a similar conductance. The 
lower dissipation is the result of the suppressed density of 
states, see Eqs. (|4.8[) . (|5.7[) and (|5.8[) and the discussion 
following these equations. 

On the other hand, we have shown that at low temper- 
atures the fluctuations of the admittance caused by the 
fluctuations of the Andreev level occupation can become 
large (see Section IVI[) . At fixed number of conducting 
channels N e , they are larger than the admittance fluctu- 
ations of a tunnel junction by a factor exp[(A — Ea)/T], 
where Ea is the energy of the lowest Andreev level, see 
Eqs. (|6.2[) and (|6.5p . In addition to that factor, already 

— 1/2 

enhancing fluctuations, their amplitude scales as N e ' . 
Josephson junctions in the existing qubit devices have 
conductance G ~ e 2 /7r. An all-metallic link replacing 
such junction would have N e ~ 1 leading to gigantic 
fluctuations. Situation is better for resonant devices de- 
signed for different applications^ where N e ~ 100, and at 
the same time the demand on the resonance frequency 
stability may be milder. 

The admittance of a single-channel junction exhibits 
a resonant behavior at frequencies uj sa 2Ea correspond- 
ing to the creation of pair of quasiparticles in the An- 
dreev level. We studied in more detail the effect of the 
AC perturbation of such frequencies on a quasiparticle 
dynamics, see Section IVIII If Ea < A/3, the system 
goes through Rabi oscillations between the empty and 
doubly-occupied Andreev level without dissipation. For 
Ea > A/3 the AC perturbation also causes excitations 
of quasiparticles from the level to the band. The dissipa- 
tion power is then non-zero and has resonant behavior, 
with the resonance width depending on the amplitude 
of the AC perturbation, Eq. (|7.8|) . We found that the 
junction inductance follows the same behavior, therefore 
the Rabi frequency can be measured in a two-tone ex- 
periment. In the case of a disordered weak link, there is 
no dissipation at the AC perturbation frequencies lower 
than 2A/3. At higher frequencies, the dissipation power 
depends non-linearly on the AC perturbation intensity, 
see Eq. (|7.11[) . Finally, it is worth noting that the pop- 
ulation of Andreev level may depend non-monotonically 
on the intensity of perturbation, see Eqs. (|7.7p and (|7.9[) . 
The population of separate Andreev levels may be stud- 
ied in experiments^ with break junctions. 



Acknowledgements 



W911NF-09-1-0369. 



Appendix A: The complete expression for ImY(oj) 



From Eq. (|2.5[) we can get the complete expression 
for the admittance, including the imaginary part. Since 
Y(uj) is analytical in the upper half-plane, ImV(u) can 
also be obtained from the expressions for ReF(w) by 
Kramcrs-Kronig relations. At zero temperature, the con- 
tributions to the Im Y(u;) corresponding to the real parts 
from Eqs. p.3[) - (|3.7p are given by 



ImYi(w) _ I 
G ttu. 



r-OO />00 

P.V. / dei I de 2 p{e 1 )p(e 2 )\z euei+e2 \ 

J A J A 



x [l-/( ei )-/(e 2 )] 



OJ — 1\ — 6 2 Ld + £i + £2 £\ + (-2 



(Al) 



Im r 2 H = y/A 2 - E\ 
G uj 



poo 

P.V. / dep(e)\z e . e+E J 

J A 



1 



f 



U) — e—E A ui + e+E A e+E A 

(A2) 



Im y 3 H (A 2 -E\){E\- A 2 cos 2 f) 



G 



ujE\ 



1 



I 



uj - 2E A uj + 2E A E A 



(A3) 



These expressions are valid for any distribution of quasi- 
particles. The case of no quasiparticles present corre- 
sponds topo = lj Pt,i,2 = and /(e) = 0. The imaginary 
parts of the last two contributions, Y^(uj) and Y 5 (uj) are 
given by: 



^M = J-p.v./ &x 

G 7TUJ 



*2p(ei)p(e 2 )|z ei , ei - e2 | 2 



x [/(c 2 )-/(ei)] 



f 



■ + - 



ImY 5 (uj) -- 
„ ,Pf+Pl 



-P2-f(e)} 



cj-ei + e 2 w + ei-e 2 ei — e 2 



P.V. / dep{e)\z^_ EA \ 

J A 



(A4) 



1 



f 



uj — e + Ea uj + e — Ea e — Ea 

(A5) 
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Appendix B: Admittance of a weak link 

Let Po(e), p-f,4,(e) and p 2 (e) be the probabilities to have 
zero, one or two quasiparticles in the Andreev level with 
energy e. The occupation factor of the continuum state 
with energy e is denoted by /(e). The admittance of a 
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disordered weak link for general occupation numbers is 

given by: a 

1 ttGA f , A 2 cos 2 | - e 2 sin 2 t. 



Y(u)=T i Y i (u) + — , (Bl) Lj | sin f | J e 2 Je 2 - A 2 cos 2 £ 

f— J OjLj A|cosf| V 2 

(B2) 

where the inductance term is: 

The real parts of the Yi terms are given by: 



ReYi(w) 
G 



RcY 2 (uj) 
G 



ui — A 



(w-2A)P.V./de 



-A 2 (l+cos0) Ve 2 -A 2 v /(w-e) 2 -A 2 ln \/ e2 / A2_C0S 2' 



A 



7T / de6(uj - e — A) 



A 2 w 2 |sinf|(w-2e) y^/A^cos 2 f ye 2 /A 2 -cos 2 f-| sin 

we- A 2 (l + cos0) v /(w-e) 2 -A 2 \/A 2— 7 



.[l-/(e)-/(w-e)], 
(B3) 



A| cos f| 

Rey 3 H 



; 2 A|sin§ |(w-2e) 



: 2 — A 2 cos 2 $ 



'wO+^-^-O], (B4) 



A 2 1 



0(w-2A|cos||)0(2A-w) — — I — --1 ^— -4cos 2 |[p (^/2)-p 2 (w/2)], (B5) 



ReKi(w) 



G 



de0(e + w-A)[/(e)-/(e + w)] 



Ve 2 - A 2 y/(e + w) 2 - A 2 
w 2 A 2 |sinf |(2e + w) 



(e+w) - A 2 (l +cos 0) ^ V /(e+w) 2 /A 2 -cos 2 | + |sin|| ^ A 2 (l+cos</.) ^ ^/e 2 /A 2 -cos 2 | + 1 sin § | 



v /(e+w) 2 /A 2 -cos 2 | v /(e+w) 2 /A 2 -cos 2 |-|sin|| ^e 2 /A 2 -cos 2 £ v /e 2 /A 2 -cos 2 § - 1 sin f 



(B6) 



G 



7r / de6»(e + w - A) 



we + A 2 (l + cos0) v /(w + £ ) 2 -A 2 VA 2 -e 2 
w 2 A|sinf|(2e + w) / £ 2 - A 2 cos 2 | 



+Pa(e)-f(u+e) 



(B7) 



From these we can also find imaginary parts using consider two opposite limits, w <C T: 
Kramers-Kronig relations. 

Assuming Boltzmann distribution /(e) = Ae~ e ^ T both RcKdw) ~ x cot 2 — \f~KjT (B8) 

below and above the gap, for frequencies less than \/2tt 2 

2A|cos(0/2)| the only contribution to ReF(w) comes _ 7r 2 G 2 w 

from dHSJ) and (|B7)l . In the limit w,T < A we get: ReF 5 (w) cot 2 ' ^ 



ReF 4 (w) « -=x gp cot -y — (1-e ' )Z7(-) , 



G 2 Tv^ w 



tt 2 G 2 ^ VTA w w ReY 4 (w) « — ^=i w cot^— 375-ln-, (BIO) 

V2^ 9P 2 w 2T n 2T> 

Rey 5 (w)«^x 9J3 cot 2 ^^=e" /T . (Bll) 



7rG 9 6 vTA , 

^ cot — — ,="/-' 

where C/(x) = / °° dte- xt y/x(l + x) ln(l + l/x) and h{x) — av ~; 2 V / 2 9P 2 w 3 / 2 



is the modified Bcsscl function of the first kind. Now 
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